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(57) Abstract 

A communications system employing the Discrete Multitone (DMT) line 
code divides the used frequency band W into a large number of subchannels, 
each with a bandwidth A/-W/W. The signal in each subchannel is independently 
modulated at the symbol rate 1/A/ Because of noise and the non-ideal 
characteristics of the transmission channel, the signal-to-noise ratio (SNR) 
varies for different subchannels across the frequency band. Each subchannel 
is used to transmit an integer number of bits in a two-dimensional signal 
constellation, the susceptibility to noise increasing with the size of the 
constellation assuming that the same signal energy is used. In order to use 
the transmission channel optimally, subchannels are assigned different numbers 
of bits depending on their specific SNR values, i.e., a subchannel with a high 
SNR value is assigned a large number of bits and vice versa. An optimal bit 
allocation is determined that ensures a certain bit error rate (BER) is obtained 
for all subchannels. In order to simplify numerical computations, thresholding 
is employed to make initial bit assignments. Further simplification of the bit 
loading procedures is obtained by ending bit allocation optimization efforts after 
a targeted BER is achieved. 
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REDUCED COMPLEXITY BIT ALLOCATION TO SUBCHANNELS IN A 
MULTI-CARRIER. HIGH SPEED DATA TRANSMISSION SYSTEM 



FIELD OF THE INVENTION 

The present invention relates to multi-carrier, high speed data 
5 transmission and reception, and more particularly, to a method and system for simply 
and cost effectively allocating optimal bandwidth to each of the multi-cam 

BACKGROUND OF THE INVENTION 

The general principle of transmitting data by dividing it into several 
interleaved bit streams, and using these bit streams to modulate several carriers is well 
jo known. Modems using this principle employ what is referred to as Discrete Multitone 
(DMT) modulation and demodulation. These types of modems are being used or 
considered for use in such applications as cellular radio and High rate Digital 
Subscriber Lines (HDSL). 

In a Discrete Multitone system, the input bit stream is first serial-to- 
15 parallel converted. The parallel output is then grouped into N groups of bits 

corresponding to the number of bits per symbol. Portions of bits are allocated to each 
DMT carrier or subchannel. The power transmitted over each subchannel corresponds 
to the number of bits transmitted over each DMT subchannel . 

Like a single carrier system, a multicarrier data transmission system tries 
20 to maximize a data transmission or bit rate for a particular Bit Error Rate (BER) and a 
particular Signal to Noise Ratio (SNR). In theory, the total bit rate (transmitted over all 
of the sub-channels) is maximized if the bit error rates for all the subchannels are equal. 
The concern is how best to allocate bits to each subchannel to attain the same BER over 
all of the subchannels. 
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Fig. 1 shows an example Discrete Multitone (DMT) communication 
system in which the present invention may be advantageously employed. Transmitter 
10 includes a serial-to-parallel converter 14, a multicarrier modulator 16, and a 
pretransmit processor 1 8. Receiver 12 includes a post channel processor 20, a 

5 multicarrier demodulator 22, and a parallel-to-serial converter 24. The transmitter and 
receiver are linked in this example by a digital subscriber line (DSL) or other form of 
communication channel 26. Serial input data at a rate of btotal^T bits per second are 
grouped by converter 1 4 into blocks of btotal bits for each multicarrier symbol, with a 
symbol period of T. The btotal bits in each multicarrier symbol are used to modulate N 

10 separate carriers in modulator 16 with bj bits modulating the r^ 1 carrier. 

A preferred example embodiment uses an Inverse Discrete Fourier 
Transform (IDFT) during modulation to generate N s time-domain samples of a transmit 
signal for each block of btotal bits, where N s is preferably equal to 2N. The : 
corresponding multicarrier demodulator performs a Discrete Fourier Transform (DFT), 

15 where bj bits are recovered from the r*h carrier. As depicted in Fig. 2, the carriers or 
subchannels in a DMT system are spaced 1/T Hz apart across N/T Hz of the frequency 
band. More detailed discussion of the principles of multicarrier transmission and 
reception in general is given by J. A. C. Bingham in "Multicarrier Modulation for Data 
Transmission: An Idea Whose Time Has Come", IEEE Communications Magazine^ 

20 Volume 28, Number 5, pp. 5-14, May 1990. 

The general structure of a Discrete Multitone Modulation DMT system is 
illustrated in Fig. 3, where {Xo^Xj,..., X]sj.i } are the original, complex, input data 
symbols, {xfc} is the modulated data sequence (before cyclic prefix), {h^j is the 
discrete-time channel response, {n^} is the additive noise sequence, {y^} is the 
25 received sequence (after the removal of cyclic prefix), and { 3c q, ?j, r N- 1 } are the 
decoded, complex data symbols. The pj's and pj*'s in Fig, 3 are known as the 
modulating and the demodulating vectors, and preferably they are chosen to be 
orthonormal. 
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Adding the cyclic prefix is a discrete-time technique used to eliminate 
interblock interference (IBI) in the DMT system. The independent modulating and 
demodulating vectors in Fig. 3 are the IDFT and the DFT vectors, given by the 
following pair of relationships: 

: 1 27V-1 -J 2mm 

I 2jV-I f2mwi 

A DMT system with N complex subchannels in the frequency domain 
requires a DFT size of 2N, and the forced conjugate symmetry in the frequency domain 
results in the desired real-valued time domain samples. In the preferred embodiment, 
10 IDFT and DFT are implemented with the well known IFFT and FFT algorithms using 
the approach disclosed in commonly-assigned "Method and Apparatus for Efficient 
Computation of Discrete Fourier Transform and Inverse Discrete Fourier Transform," 
U.S. patent application Serial No. 08/887,467, (attorney docket 14 10-340), filed on July 
2,1997. 

15 Fig, 4 is a more detailed block diagram showing the principal operative 

components of a basic DMT transmitter 30 and a basic DMT receiver 32 connected 
through a channel 34. Serial input data are grouped into blocks, converted to a parallel 
form, and appropriately encoded by an encoder 36. Each one of the N active 
subchannels contains a number of bfc bits allocated to that subchannel. The total 

20 number of bits in a DMT symbol is therefore: 

n \ - t 

B to<al = Y. b k (3) 
k=\ 



(i) 

(2) 
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Each of the N bit groups is mapped by the symbol into a two-dimensional 
signal constellation. The relation between the number of bits in each block of bits bfc 
and the necessary number of signal points in a constellation M is given by: 

M k =2 b > (4) 

5 The signal points in a constellation are typically arranged in a rectangular 

pattern for even values of bfc and a cross for odd values of bfa Provided that the same 
signal energy is used for all constellations, the distance between neighboring signal 
points is reduced with increasing constellation size resulting in increased BER with the 
same signal-to-noise ration (SNR). 

10 The output from the encoder is N complex numbers, one for each group of 

bits, which are then fed to a device that calculates the Inverse Discrete Fourier 
Transform (IDFT). The output is a real sequence that can be considered the 
superposing of TV modulated orthogonal carriers spaced A/ apart. 

The parallel IFFT output is converted back to a serial data stream by a 
15 converter 40. The digital modulated data stream is cyclically prefixed, converted to 
analog form by a digital-to-analog converter (DAC) 42, low-pass filtered at 44, and 
passed through a D.C. isolating transformer 46 during pre-transmit processing to 
produce an analog signal transmitted over a transmission channel 34. 

At the receiver end, the received analog signal is passed through a D.C. 

20 isolating transformer and low-pass filter 48, converted to digital form by an analog-to- 
digital converter (ADC) 50, time domain pre-equalized by a finite impulse response 
(FIR) filter 52 to limit the effective memory of the channel, and stripped of the cyclic 
prefix during post-receive processing in converter 54. The resulting digital signals are 
demodulated by 2N FFT operation 56 and converted to parallel frequency domain 

25 signals. Since the amplitude vs. frequency and the delay vs. frequency responses of the 
channel are not necessarily constant across the entire used frequency band, the received 
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signal will differ from the transmitted signal, and the parallel inputs to the decoder 58 
will differ from those parallel outputs from the encoder 36 ? One form of equalization 
used to compensate these differences is a frequency domain equalizer (FEQ) which 
individually adjusts for the attenuation and delay of each of the carriers immediately 
5 before the parallel frequency domain signals are passed to the decoder. The frequency 
domain equalized signals are appropriately decoded and converted back to a serial form 
by the decoder 58. Ideally, the detected output serial data from the decoder 58 is the 
same as the input serial data to the encoder 36. 

In a conventional single-carrier digital transmission system, such as a 
10 Quadrature Amplitude Modulation (QAM) system implemented with an equalizer or a 
encoder, the transmission bandwidth is determined by the symbol rate and the carrier 
frequency of the system. Unfortunately, variable symbol rate, single-carrier systems are 
still impractical to implement from a complexity standpoint with today's technology, 
and even if they are implemented, the granularity of possible symbol rates are typically 
15 very coarse. This is a direct consequence of the fact that for a fixed data rate, the 
symbol rate can only change in discrete multiples of b S y m b 0 ]/(b S y m b 0 l±l) 5 where 
b S y m bol i s ^ e number of bits transmitted by each data symbol, provided that only 
signal constellations with integer numbers of bits are used. 

Multicarrier modulation, however, offers much more flexibility in 
20 granularity because it acts on block symbols that consist of a much larger number of 
bits per symbol over a large number of carriers. Since different numbers of bits can be 
transmitted through the different carriers (subchannels), the multicarrier transceiver can 
control of the transmission bandwidth. As a result, more data bits are transmitted 
through the better carriers, i.e., having a greater Signal to Noise Ratio (SNR), and less 
25 data bits are transmitted through the worse carriers, i.e., having a lower SNR as shown 
in Fig. 5. In this manner, optimal performance can be achieved, in terms of either 
maximizing total data throughput for a fixed system performance margin or maximizing 
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overall system performance margin for a fixed target data rate under a fixed, total input 
power constraint and target bit-error-rate. 

Various bandwidth/bit allocation optimization techniques for multicarrier 
modulation have been proposed for different applications. The theoretical optimal 

5 energy allocation is the well known water-pouring energy distribution as described for 
example in J.G. Proakis, (1995), Digital Communications 3 rc * edition, McGraw-Hill, 
New York. Conceptually, the water pouring energy distribution model is obtained by , 
first inverting the channel SNR curve shown in Fig. 5 in the frequency domain and then 
pouring the available energy into the resulting curve. The final "water level" at each 

io particular frequency /carrier/subchannel is the optimal amount of energy to use at that 
particular frequency/carrier/subchannel. The number of bits transmitted by each carrier 
is determined by the maximum number of bits that can be supported by each carrier, 
given the water-pouring energy distribution and the desired bit-error-rate of each 
carrier. While the water-pouring technique yields the optimal bit and energy 

15 allocations, and therefore the optimal transmission bandwidth, it is typically complex 
and impractical to implement. 

Another multicarrier bandwidth optimization technique designed for 
voice-band modem applications, where only constellations with an integer number of 
bits are used, was proposed by D. Hughes-Hartogs and described in "Ensemble Modem 

20 Structure for Imperfect Transmission Media", European Patent Number 0 224 556 Bl 
(published October 1991). The basic idea behind Hujghes-Hartogs' algorithm is to 
incrementally add, one bit at a time, the amount of data to be transmitted in each 
. multicarrier symbol until a desired data rate or power constraint is satisfied. 
Furthermore, this algorithm always chooses the subchannel that requires the least 

25 amount of incremental power to support an additional bit at the given desired bit-error- 
rate. Hughes Hartogs' algorithm, however, is typically too slow in convergence, 
making it especially unsuitable for applications where there are large numbers of 
carriers and bits to be supported by each multicarrier symbol as pointed out in P.S. 
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Chow, J.M. Cioffi, J.A.C. Bingham, "A Practical Discrete Multitone Transceiver 
Loading Algorithm for Data Transmission Over Spectrally Shaped Channels/' IEEE 
Transactions on Communications, Vol. 43, No. 2/3/4, pp. 773 775, 1995. Also, the 
channel response is assumed to be fixed, i.e., time-invariant after training, which is not 
5 a good assumption for many applications, e.g., mobile radio communications. 

In order to achieve faster convergence at the expense of a slight 
performance degradation, Chow et al. proposed a three-step bit loading algorithm. 
First, an optimal system performance margin is found; second, convergence is 
guaranteed with a suboptimal loop; and third, the energy distribution is adjusted for 
10 each subchannel if the transmission system is constrained by the total transmit power. 

However, Chow's bitloading algorithm requires computationally complex 
steps best run on a microprocessor and a number of iterative elements that may be 
excluded. Both requirements increase the system cost and reduce speed. The bit 
loading procedure in accordance with the present invention minimizes computational 
15 complexity and permits execution by specific hardware which is faster than 

programmed microprocessor driver solutions. Also those parts of the algorithm that are 
best run on a microprocessor are of such basic character that only a rudimentary, 
inexpensive device is necessary. 

It is an object of this invention therefore to provide optimal bitloading 
20 using less computational power and at a fast rate. 

SUMMARY OF THE INVENTION 

Bits are allocated among subchannels in a multi-carrier modulation 
environment, under the constraint of a targeted bit error rate (BER) as well as a targeted 
number of bits per symbol. Thresholding is used to make initial bit assignments to 
25 subchannels in order to simplify numerical calculations. If the targeted BER cannot be 
achieved with the targeted number of bits, the overall channel is insufficient, and a 
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higher BER or a lower transmission rate must be considered. To further increase the 
efficiency of the process, a "block-reduction" step may be employed which reduces by 
one the number of bits allocated to each active subchannel if the total number of 
allocated bits exceeds the number of targeted bits by more than the number of active 
5 subchannels. Further fine tuning may also be performed. 

The simplified bit allocation procedures in accordance with the present 
invention may be performed quite rapidly to permit adaptative allocation of bits in real 
time to accommodate, for example, changing channel conditions. In addition, the 
simplified procedures permit implementation using dedicated hardware for improved 
10 speed. 

More specifically, in a multi-carrier communications system .with a 
transmitter and a receiver linked together by a communications channel affected by 
noise and other channel impairments, data bits to be transmitted over the channel at a 
prescribed data transmission rate and bit error rate are allocated to plural subchannels. 

15 A series of thresholds is calculated for each subchannel with each threshold 

corresponding to a preset number of bits to be allocated to that subchannel. Each 
subchannel is associated with a signal constellation where bits allocated to the 
subchannel are mapped to signal points in the signal constellation. The subchannel 
thresholds are calculated based on the prescribed error rate, an average energy and size 

20 of the subchannel signal constellation. The threshold calculations need only be 
performed once, preferably before the transmitter/receiver is brought on-line. 

Subsequently, after the transceiver is brought into operation, an error 
estimate such as a noise standard deviation is obtained for each subchannel. That error 
estimate is compared to the thresholds calculated for the subchannel. The number of 
25 bits is then allocated to each subchannel corresponding to one of the thresholds 
exceeded by the obtained error measurement. The error estimate comparing and 
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obtaining procedures as well as the bit allocation may be performed in real time 
whenever the channel condition changes. 

A determination is made whether a total number of bits allocated to all the 
subchannels is greater than or equal to a total target number of bits associated with the 

5 prescribed data transmission rate. If not, a decision is made whether to reduce the 
prescribed data transmission rate or to reduce the prescribed bit error rate. Otherwise, 
the number of bits allocated to one or more of the subchannels may be reduced until the 
total number of allocated bits is nearly equal to the total target number of bits. In - ' 
particular, the number of bits to be allocated to each active subchannel may be reduced 

10 by one until the excess number of bits is less than the number of active subchannels. 

When a remainder of the excess number of bits is less than the number of 
active subchannels, the remainder may be further reduced by reducing the number of 
bits allocated to selected ones of the active subchannels. Bits are subtracted based on a 
magnitude of a difference between the error estimate and the exceeded threshold for 
15 each active subchannel. Those active subchannels having larger difference magnitudes 
are selected first for the bit reduction followed by subchannels having smaller 
difference magnitudes. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 is a block diagram illustrating a basic prior art multicarrier 
20 transmitter, communication channel and receiver; 

Fig. 2 is a diagram conceptually illustrating frequency domain discrete 
multitone modulation (DMT); 

Fig. 3 is a diagram schematically depicting the basic structure of a 
discrete multitone modulation system; 
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Fig. 4 is a block diagram depicting the principal components of a discrete 
multitone (DMT) modulation communication system; 

Fig. 5 illustrates DMT bit allocation in the frequency domain; 

Fig. 6 illustrates two-dimensional 4 and 16 QAM signal constellations; 

Fig. 7 shows one example configuration implementing the thresholding 
operation and bit allocation procedures in accordance with one example embodiment of 
the present invention; 

Fig. 8 is a diagram providing a specific example threshold and bit 
allocation operation, based on the example configuration shown in Fig. 7; 

Fig. 9 shows an alternative example configuration for implementing the 
thresholding and bit allocation operations in accordance with the present invention; and 

Fig. 10 illustrates bit loading procedures in flowchart form in accordance 
with the present invention. 

DETAILED DESCRIPTION OF THE INVENTION 

In the following description, for purposes of explanation and not 
limitation, specific details are set forth, such as particular hardware, applications, 
techniques, etc. in order to provide a thorough understanding of the present invention. 
However, it will be apparent to one skilled in the art that the present invention may be 
practiced in other embodiments that depart from these specific details. In other 
instances, detailed descriptions of well known methods, protocols, devices, and circuits 
are omitted so as not to obscure the description of the present invention with 
unnecessary detail. 

The example apparatus that may be used in implementing the present 
invention may include the transmitter and receiver represented by the block diagram of 
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FIG. 4 already described above that in addition includes hardware and/or software for 
implementing the novel bit allocation procedures according to the present invention. 
The theoretical foundation for those bit allocation procedures now follows. 

Theoretical Foundation 

5 The bitloading procedures in accordance with the invention allocates bits 

to be transmitted in such a way that each subchannel satisfies the demands on 
probability of a symbol error, Pj^, within a system performance margin. This bit 
allocation is based on a noise variance, cr , in each subchannel and prior knowledge of 

the average energy in each subchannel signal constellation, £ av , which together 
10 effectively provide the signal-to-noise ratio SNR for each subchannel. Any additional 
safeguard margin of bits is distributed evenly among the subchannels. 

According to J. G. Proakis, (1995), Digital Communications 3 r d edition, 
McGraw-Hill, New York, the symbol error probability Pfrj for a rectangular 
constellation in Additive White Gaussian Noise (A WGN) is tightly upper bounded by: 

where a number of constellation points M per subchannel is defined as 
A/ = 2 V ; 6 = l,2,....,6 max is the number of bits conveyed per each signal point in the 

constellation; £ av is the average energy of the subchannel constellation and is a known 
design parameter set for each subchannel constellation; and N 0 is the single-sided 
20 power spectral density of the noise for a received signal constellation point (at the 
receiver) and is defined as N Q = 2cr 2 \ where cr is the received noise variance per 
real/imaginary dimension. The noise variance can be visualized (in a simplified 
manner) as the variance of the symbol error, i.e., between the signal point received and 
the closest constellation signal point in the complex plane. 
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The Q-fimction in equation (5) relates the signal-to-noise ratio (roughly 
defined as e av I N 0 ) to the symbol error for a particular constellation size and is 
defined mathematically as: 

= -^-]'e^du (6) 

5 where u is the integral variable. In other words, the Q-function of a certain value x is 
calculated by solving the integral between the limits x and infinity. 

Even though Yy[ and BER are not identical, this analysis assumes that the 
symbol error probability and bit error rate (BER) are approximately the same because 
one symbol usually corresponds to one or more bits. Certain symbol encoding 
10 techniques, (i.e., Gray coding where adjacent symbols containing the most likely errors 
only differ by one bit), may be employed to ensure that the bit error probability and the 
symbol error probability are very similar. 

In most practical applications, the bit error rate (and in this application the 
symbol error rate Pj^) is a known system design constraint. To achieve the same BER 
15 (or Pm) for all subchannel signal constellations and thereby optimally allocate bits to 
each subchannel under a constraint that the average energy of all subchannel 
constellations is the same, the denominator of equation (5), (M-l) No must be 
maintained constant, i.e., 

(M~l)iV 0 = 2^(2 6 ~l) = C (7) 

20 The constant C depends on the specific BER to be achieved and may be calculated 
according to: 

(8) 



C = 3f„ 



Q 



-1 
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Solving equation (7) for the number of bits to be allocated, b, provides the 
number of bits that may be assigned to a subchannel given a certain BER and the noise 
variance c\ on that particular subchannel, i.e., the bit allocation for each subchannel k 
corresponding to b k : 

where a\ is the noise variance per real (I) and imaginary (Q) dimension on 
subchannel k. In a preferred implementation, the bit loading values are limited to 
integer constellations, e.g., by rounding them to integer constellation sizes which affects 
the total number of allocated bits. 

10 In equation (7), both o and b may vary. The subchannel noise variance 

cr k varies because of changing noise on the subchannel. Therefore, the number of bits 
. 6^ allocated to that subchannel must be varied to compensate for the subchannel noise 
changes to ensure that equation (7) is satisfied as a constant and equal to C. If the 
subchannel noise variance o\ increases, b^ must decrease, and fewer bits are sent over 

15 the subchannel. If & k 2 decreases, b^ can be increased with more bits being sent. 

As already described, for a particular system design, the parameters Pm 
or BER, Q, and s av are preset/known. Therefore, only the noise variance a 1 needs to 
be determined for each subchannel to solve for the optimal bit allocation per subchannel 
according to equation (9). To obtain the subchannel noise variance a\, the noise 
20 standard deviation o k is determined based on a series; of error measurements for 

symbols received on the subchannel. The subchannel noise standard deviation a k is the 
input to the bit loading procedures employed in the present invention and corresponds 
to the mean of the received symbol error magnitude (attributed to noise) for each 
subchannel. 
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The following described how the noise standard deviation estimate a* 
may be determined for eaeh subchannel k. The received symbol error is the vector 
difference between a received signal point and the closest signal constellation point for 
a particular subchannel. Thus, the received symbol error is a complex number. To 

5 estimate the deviation of the error, either the real (I) or the imaginary (Q) part of the 
received symbol error may be estimated for a theoretical ideal channel with truly 
Additive White Gaussian Noise (AWGN) because the I and Q components would have 
equal variance, i.e., the actually received signal points would be equally distributed in a 
circular pattern around the signal constellation point. In reality, the noise variance may 

io be different in the I and Q dimensions, i.e., riot a circular pattern, resulting in poor noise 
deviation estimates if noise in only one of the two dimensions is measured. 

Thus, it is preferred that the noise magnitudes in both the I and Q 
dimensions be estimated. The received symbol error magnitude can be expressed as: 



R = + x \ (10) 

15 where X { and X 2 are zero-mean, statistically-independent, gaussian-random variables, 
each having a variance a 2 corresponding to the noise components in the complex I- and 
Q-dimensions. According to J.G. Proakis, ( 1995), Digital Communications 3 r d edition, 
McGraw-Hill, New York, the random variable R is Rayleigh distributed with an 
expected (E) mean value: 

20 a = E(R) = ^a (11) 

From equation (11), the measured noise standard deviation estimate a k for each 
subchannel is preferably corrected using E(R), (i.e., by dividing the measured a k by 



y , in order to obtain the true value of a k . 
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Example Embodiment of Subchannel Bit Allocation Procedures 

The subchannel bit allocation procedures in accordance with the present 
invention take as an input vector the corrected noise standard deviation estimate a k for 
each subchannel k. As described above, the main goal in optimally allocating bits is to 

5 maintain a constant, predetermined symbol error probability P]yj/BER across each of the 
subchannels. Towards that goal, the invention uses thresholding to make initial bit 
assignments rather than using conventional complicated numerical computations. 
Further simplification is achieved by ending bit allocation optimization efforts after a 
target Pjyj/BER is obtained. In another example embodiment, the preferred bit 

io allocation procedures includes three basic steps that each converge within a finite 
amount of time. . . ? , 

Thresholding 

A corrected, standard deviation of the noise a k for each subchannel k is 
estimated using a number of consecutively received symbols as described above in 

15 conjuntion with equations (10) and (1 1). A maximum number of bits to be allocated for 
each subchannel is then determined by comparing the estimated noise standard 
deviation a k measured for each subchannel A: to a series of thresholds r previously 
calculated and stored for each subchannel. Thus, rather than having to perform the 
complicated calculations in equations (8) and (9) including logarithms, exponentials, 

20 multiplications, etc., the present invention simply compares the estimated noise value 
measured for each subchannel to a number of thresholds. 

The thresholds themselves are calculated only once, (e.g., off-line at a 
factory), and are subsequently used to establish a correspondence between a subchannel 
noise standard deviation value a k and a number of allocated bits bfc for that subchannel. 
25 In other words, when the subchannel noise standard deviation value exceeds a certain 
limit, the number of bits that allocated to that subchannel changes. The subchannel 
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thresholds % are obtained by solving equation (7) for a and inserting equation (8) 
which yields: 



3*. 



Q 



4 



2(2 9 - l) 



(12) 



where g signifies the discrete values at which the subchannel signal constellation size 
changes. In effect, the subchannel constellation size changes with the measured 
subchannel noise to maintain a constant BER. The thresholds define the values of 
measured subchannel noise standard deviation a ' where the constellation size changes 

as given by: 

.13 1 

where b^ is the maximum number of bits to be transmitted per subchannel, e.g., a 
maximum number of 12 bits. Once the thresholds are determined (preferably off-line) 
using equations (12) and (13) and the known values e mfi and Pm and the function Q(x) 
for each subchannel, they are stored in memory along with a corresponding number of 
bits associated with each threshold. 

If the measured subchannel noise is so strong that not even the smallest 
constellation may be transmitted with the required bit error rate, that channel is not 
used, resulting in a zero bit allocation. Any subchannel with a bit allocation of zero 
reduces the number of active subchannels by one. 

The thresholding process is very simple compared to the complex 
calculations required by conventional bit loading algorithms. The thresholds are 
calculated following equations (12) and (13) and stored in memory. The received 
symbol mean error magnitude is measured for each subchannel and used as an estimate 
of the subchannel noise standard deviation a k . The subchannel bit allocation b^is 
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associated to an interval between the two thresholds in which the value of the noise 
standard deviation estimate <j k falls. The thresholding/comparison may be 
implemented in software using a rudimentary data processor, or preferably in relatively 
simply hardware to increase operational speed. 

5 Fig. 7 is block diagram illustrating the thresholding operation. The 

measured subchannel noise standard deviation a k is input to a bank of parallel 
comparators 70. The other input of each comparator 70 is one of a number of 
subchannel thresholds r calculated in accordance with equation (12) for different 
values of q as indicated in equation (13). In an example embodiment where b max is 12, 

io there are twelve comparators 70 arid twelve different stored threshold values r, - t n 72. 
The combined outputs of each comparator 70 for a "bit vector" that is input to an 
encoder 74. The least significant bit (LSB) of the bit vector is fixed to zero. In each 
comparator 70, if the input noise deviation a k is less than the specific threshold value 
r, the comparator 70 generates a logical u i" at its output. Alternatively, if the input 

15 value is greater than the threshold, the comparator 70 generates a logical '."0" at its 
output. 

As illustrated in Fig. 8, the least significant bit in the comparator input bit 
vector supplied to the encoder 74 is in the zero position, and the encoder 74 outputs a 
zero if the bit vector contains only zeros. In effect, the specific comparator whose input 
20 threshold is first exceeded by the input noise standard deviation r k determines the 
number of bits indicated by encoder 74 to be allocated to that subchannel. In the case 
shown in Fig. 8, the smaller the noise standard deviation value, the larger the number of 
bits that may be allocated to the subchannel. 

Of course, the thresholding operation illustrated in Figs. 7 and 8 is only an 
25 example. Another example implementation is shown in Fig. 9 where a single 
comparison device 84 is used. The measured noise standard deviation r k for a 
subchannel is stored in an input register 80 and input to a comparator 84. Digital 
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threshold values r calculated for this subchannel are stored in threshold, memory 82 
and are sequentially input to the digital comparator 84. The resulting sequential 
comparator outputs are passed through a demultiplexer in order to write different 
positions in the bit vector register 88. The encoder 90 analyzes the bit vector stored in 
5 register 88, and outputs the corresponding number of bits to be allocated for that 
subchannel based upon the thresholding comparison operation. 

: If the total number of bits allocated to all subchannels in the thresholding 
step is less than the target number of bits, the desired bit rate cannot be achieved while 
maintaining the prescribed bit error rate because of increased noise or other factors. 

io Consequently, the data transmission rate or the BER performance must be reduced. A 
decrease in transmission rate means allocating ai fewer number of bits to each 
subchannel. Alternatively, another set of thresholds t corresponding to a higher BER 
may be employed. Thus, several sets of threshold values for different BERs may be 
calculated using equations (12) and (13) and stored in memory. However, the former 

15 decrease in transmission rate/allocated bits is preferred. 

Block Reduction 

If the total number of bits allocated across all subchannels, exceeds the 
target number of bits, excess bits may be allocated to the active subchannels such that 
the data transmission rate can exceed the required data rate and still fulfill the 
20 prescribed bit error rate. If that goal is desired, the bit allocation procedures may end 
here. However, further bit allocation refinements may also be desirable. For example, 
the transmitted bit rate may be reduced to the prescribed bit rate to achieve certainty as 
well as an extra immunity to noise. 

I f such further refinement is desired, a decision is made if the total 
25 number of allocated bits exceeds the target number of bits by more than the current 
number of active subchannels, i.e., subchannels that have not been set to zero. If so, 
one bit is subtracted from the initially allocated number of bits to be transmitted on each 
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active subchannel. Of course, any subchannels that have zero bit allocations reduces 
the number of used subchannels by one. 

This "block reduction" step is repeated until the number of allocated bits 
in excess of the overall target number of bits is less than the number of active 
5 subchannels k. The reduced number of allocated bits corresponds to a lower overall 
transmission rate with the margin evenly distributed among the subchannels. 

Fine Correction 

A further fine correction step may be employed to compensate for the 
difference between the number of bits remaining to be allocated and the target number 
10 of bits after the block reduction step is completed. For each active subchannel, the 
relation between the estimated noise standard deviation a k and the threshold of its 
subchannel constellation is calculated as: 

where thr[b k ] signifies the operation of comparing the noise standard deviation a k of 
15 subchannel k to the precalculated thresholds for that subchannel. The distance 
represents the difference between the noise standard deviation value a k and the 
associated threshold r where the bitloading value changes. The subchannel having the 
largest distance is the first subchannel to have a bit removed from the groups of 
allocated bits. 

20 Because bits are subtracted from subchannels with the largest distance d k 

the subchannel most sensitive to noise in its present constellation is advantageously 
reduced first. If the bit allocation of a subchannel that is not the most sensitive to noise 
is reduced, its immunity to noise is boosted. Unfortunately, the most sensitive 
subchannel remains unchanged thereby becoming the "weakest link". The distance 
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calculation ensures that the weakest links are the first to have bits removed then a 
relatively strong subchannel. 

This task may take more or less time consuming depending on the total 
number of subchannels to be reduced and on the choice of sorting algorithm. For 
example, assume there are 1024 subchannels and that after the block reduction step 
there are an excess of 1023 bits allocated, i.e., 1023 bits more than the targeted number. 
The task then is to reduce 1023 in order of "largest distance first" by finding the 1023 
largest elements in the distance vector. If, on the other hand, there is an excess of only 
10 bits, only the 10 largest elements in the distance vector need be found which takes 
less time. The sorting algorithm which sorts the distance vector elements from largest 
to smallest also affects this process. A fast sorting algorithm may be more complex to 
implement, require a sophisticated microprocessor, etc. 

A flow chart diagram outlining example procedures for implementing the 
thresholding, block reduction, and fine correction steps is now described in conjunction 
with the bit allocation routine (block 100) illustrated in Fig. 10. An initial series of 
operations is performed usually at the time the device incorporating the present 
invention, such as a modem, is being manufactured. More specifically, a series of 
threshold values r(q) is calculated for each subchannel in accordance with equations 12 
and 13 (block 102). In addition, each threshold value is correlated to a preset number of 
bits to be allocated for that subchannel. The threshold values and correlated numbers of 
bits are then stored. As mentioned above, different sets of threshold values 
corresponding to bit allocations may be calculated for several different bit/symbol error 
rates. The calculations and procedures performed in block 1 02 are generally performed 
only once. 

In terms of determining the actual number of bits to be allocated for each 
subchannel, the noise deviation estimate a k is measured for each subchannel and 
corrected in accordance with equation (11) (block 104). In a preferred embodiment, the 
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actual correction factor is incorporated in the one time threshold calculation process in 
block 102 since the correction factor is known. In addition, the number of active 
subchannels N sc is set to all of the subchannels (block 106). 

The noise standard deviation estimate for each subchannel a k is 
5 compared to the thresholds r predetermined for each subchannel from smallest to 
largest magnitude. When the noise standard deviation first exceeds one of the 
thresholds;, the stored number of bits bfc corresponding to that first exceeded threshold is 
allocated to that subchannel (block 108). For any subchannel where no bits can be 
allocated, the subchannel is deemed to be inactive, and the number of active 
io subchannels N sc is decremented. 

A decision is made whether the sum of the bits allocated to all of the 
subchannels is greater than or equal to a target number of bits B tar g et (block 110). If 
not, the current channel conditions will not provide the targeted data bit rate at the 
targeted bit error rate. A decision must be made then to decide whether to change the 
15 bit error rate requirements or the bit transmission rate requirements (block 1 12). 
Depending upon whether such reduced service requirements are acceptable, control 
may return to block 102 or block 104. 

On the other hand, if the total number of bits allocated across all 
subchannels is greater than the target number of bits, another decision is made whether 

20 the number of bits that can be allocated to the subchannels in excess of the target 

number of bits exceeds or is equal to the number of active subchannels (block 1 14). If 
so, one bit may be subtracted from the number of bits currently allocated to each active 
subchannel (block 1 16). Any subchannel whose group of allocated bits is reduced to 
zero after the subtraction operation causes a decrementing of the number of active 

25 subchannels N sc . This is an iterative process with control returning to decision 
block 1 14 via the shown feedback path. 
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When the number of allocated bits in excess of the target number is less 
than the number of active subchannels N sc , the distance for each subchannel k is 
determined in accordance with equation (14) (block 1 1 8). In order of decreasing 
distance values for each subchannel k, one bit is subtracted from each active subchannel 
5 until the sum of the allocated bits equals the target number of bits (block 120). 

While the present method may be performed under relatively static 
channel conditions, it may also be advantageously employed in dynamic channel 
conditions. For example, the noise on a particular subchannel may change over time. 
Accordingly, an inquiry is made in decision block 122 whether a subchannel condition 
TO has changed. If not, continued periodic inquiries may be made or the method may 
simply end. However, if the subchannel has changed, control returns to block 104 to 
obtain the updated noise deviation estimate for one or more subchannels affected by the 
condition change, and the method repeats. 

The present invention may, of course, be applied to any 
15 telecommunications application that employs multi-carrier type modulation, e.g., DMT, 
etc. Thus, while the invention has been described in connection with what is presently 
considered to be the most practical and preferred embodiment, it is to be understood 
that the invention is not to be limited to the disclosed embodiment, but on the contrary, 
is intended to cover various modifications and equivalent arrangements included within 
20 the spirit and scope of the appended claims. 
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WHAT IS CLAIMED IS: 

1. In a multi-carrier communications system with a transmitter and a 
receiver linked together by a communications channel affected by noise and other 
varying channel impairments where data bits to be transmitted over the channel at a 

5 prescribed data transmission rate and bit error rate are allocated to plural subchannels, a 
method comprising the steps of: 

calculating for each subchannel a series of thresholds corresponding to a preset 
number of bits to be allocated to the subchannel; 

obtaining an error estimate on each subchannel; 
io comparing the error estimate for each subchannel to the thresholds calculated for 

the subchannel; and 

allocating a number of bits to each subchannel corresponding to one of the 
thresholds exceeded by the error estimate. 

2. The method in claim 1. wherein the calculating step is performed only 
15 once and the obtaining, comparing, and allocating steps are performed for an initial 

condition of the channel and again when the initial condition of the channel changes. 

3. The method in claim 2, further comprising: 

performing the obtaining, comparing, and allocating steps in real time whenever 
the channel condition changes. 

20 4. The method in claim 1, wherein the error estimate is a noise standard 

deviation measured for each subchannel. 

5. The method in claim 4, wherein the noise standard deviation is corrected 
for nonuniform characteristics. , 

6. The method in claim 1, wherein each subchannel is associated with a 

25 signal constellation where bits allocated to the subchannel are mapped to signal points 
in the signal constellation, and wherein the subchannel thresholds are calculated based 
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on the prescribed bit error rate, an average energy of the subchannel signal 
constellation, and the signal constellation size. 

7. The method in claim 1, further comprising: 

determining whether a total number of bits allocated to all of the subchannels is 
5 greater than or equal to a total target number of bits associated with the prescribed data 
transmission rate. 

8. The method in claim 7, further comprising: 

if the total number of allocated bits is less than the total target number of bits, 
reducing the prescribed data transmission rate or reducing the prescribed bit error rate. 

10 9. The method in claim 7, further comprising: 

if the total number of allocated bits is greater than or equal to the total target 
number of bits, reducing the number of bits allocated to one or more of the subchannels 
until the total number of allocated bits is nearly equal to the total target number of bits. 

10. In a multi-carrier communications system with a transmitter and a 
15 receiver linked together by a communications channel affected by noise and other 
varying channel impairments where data bits to be transmitted oyer the channel at a 
prescribed data transmission rate and bit error rate are allocated to plural subchannels, a 
method comprising the steps of: 

obtaining an error estimate on each subchannel; 
20 comparing the error estimate for each subchannel to a series of thresholds 

corresponding to a preset number of bits to be allocated to the subchannel; 

allocating a number of bits to each subchannel corresponding to one of the 
thresholds exceeded by the error estimate; and 

when a total number of bits allocated to all of the subchannels is in excess of a 
25 total target number of bits associated with the prescribed data transmission rate, 
reducing the number of bits to be allocated to each subchannel. 
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11. The method in claim 10, further comprising: 

determining an active number of subchannels having at least one allocated bit; 

and 

if the excess number of bits exceeds the active number of subchannels, reducing 
5 the number of bits to be allocated to each active subchannel by one bit until the excess 
number of bits is less than the number of active subchannels. 

12. The method in claim 11, further comprising: 

when a remainder of the excess number of bits is less than the number of active 
subchannels, further reducing the reminder by reducing the number of bits allocated to a 
10 selected one of the active subchannels. 

13. The method in claim 1 2, wherein the one active subchannel is selected 
based on a magnitude of a difference between the error estimate and the exceeded 
threshold for each active subchannel, the active subchannel having the largest difference 
magnitude being selected as the one active subchannel. 

15 14. The method in claim 13, the further reducing step further comprising: 

subtracting one bit from each active subchannel in order of decreasing magnitude 
difference. 

15. The method in claim 10, repeating the steps in claim 10 when a condition 
of the channel changes. 

20 16. The method in claim 15, wherein the steps in claim 10 are repeated in real 

time. 

17. The method in claim 10, wherein the error estimate is a noise standard 
deviation measured for each subchannel. 

18. The method in claim 17, wherein the noise standard deviation is corrected 
25 for nonuniform characteristics. 
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1 9. The method in claim 1 0, wherein each subchannel is associated with a 
signal constellation where bits allocated to the subchannel are encoded to signal points 
in the signal constellation, and wherein the subchannel thresholds are predetermined 
based on the prescribed bit error rate v an average energy of the subchannel signal 

5 constellation, and the signal constellation size. 

20. A multi-carrier communications system comprising: 

a transmitter and a receiver linked together by a communications channel 
affected by noise and other varying channel impairments where data bits to be 
transmitted over the channel at a prescribed data transmission rate and bit error rate are 
10 allocated to plural subchannels; and 

electronic circuitry performing the following tasks: 

comparing the error estimate for each subchannel to a series of thresholds 
corresponding to a preset number of bits to be allocated to the subchannel, and 
allocating a number of bits to each subchannel corresponding to one of the 
15 thresholds exceeded by the error estimate. 

21. The multi-carrier communication system according to claim 20, wherein 
the thresholds are calculated only once and the comparing and allocating tasks are 
performed when the channel condition changes. 

22. The multi-carrier communication system according to claim 20, wherein 
20 the error estimate is a noise standard deviation measured for each subchannel. 

23. The multi-carrier communication system according to claim 22, wherein 
the noise standard deviation is corrected for nonuniform characteristics. 

24. The multi-carrier communication system according to claim 20, wherein 
each subchannel is associated with a signal constellation where bits allocated to the 

25 subchannel are encoded to signal points in the signal constellation, and wherein the 
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subchannel thresholds are calculated based on the prescribed bit error rate, an average 
energy of the subchannel signal constellation, and the signal constellation size. 

25. The multi-carrier communication system according to claim 20, wherein 
the electronic circuitry determines whether a total number of bits allocated to all of the 

5 subchannels is greater than or equal to a total target number of bits associated with the 
prescribed data transmission rate. 

26. The multi-carrier communication system according to claim 25, wherein 
if the total number of al located bits is less than the total target number of bits, the 
electronic circuitry reduces the prescribed data transmission rate or reduces the 

10 prescribed bit error rate. 

27. The multi-carrier communication system according to claim 25, wherein 
if the total number of allocated bits is greater than or equal to the total target number of 
bits, the electronic circuitry reduces the number of bits allocated to one or more of the 
subchannels until the total number of allocated bits is nearly equal to the total target 

15 number of bits. 

28. The multi-carrier communication system according to claim 27, where the 
electronic circuitry performs the following further tasks: 

determining an active number of subchannels having at least one allocated bit; 

and 

20 if the total number of allocated bits exceeds the total target number of bits by the 

active number of subchannels or more, reducing the number of bits to be allocated to 
each active subchannel by one bit until the excess number of bits is less than the number 
of active subchannels. 

29. The multi-carrier communication system according to claim 28, wherein 
25 when a remainder of the excess number of bits is less than the number of active 
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subchannels, the electronic circuitry further reducing the remainder by reducing the 
number of bits allocated to a selected one of the active subchannels. 

30. The multi-carrier communication system according to claim 29, wherein 
the one active subchannel is selected based on a magnitude of a difference between the 
error estimate and the exceeded threshold for each active subchannel the active 
subchannel having the largest difference magnitude being selected as the one active 
subchannel. 

31. The multi-carrier communication system according to claim 14, wherein 
the electronic circuitry subtracts one bit from each active subchannel in order of 
decreasing magnitude difference, 

32. The multi-carrier communication system according to claim 20, wherein 
the electronic circuitry implements the tasks using a programmed data processor. 

33. The multi-carrier communication system according to claim 20, wherein 
the electronic circuitry implements the tasks using discrete hardware elements. 
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Fig 8 



THRESHOLDS 



INPUT TO 
ENCODER 



OUTPUT FROM 
ENCODER 



0.4375 

0.2082 

0.1394 
0.0876 
0.0605 

0.0423 
0.0297 
0.0167 
0.0148 

0.0104 
0.0074 

0.0052 



OBITS 



3 BITS 



4 BITS 



5 BITS 



6 BITS 



7BrTS 



8 BITS 



11 Brrs 



12Bn*S 



"0000000000000" 

I I 
MSB LSB 



1 BITS "000000000001 0" 



2 BITS "000000000011 0" 



"000000000111 0" 
"000000001111 0" 

"000000011111 0" 

"000000111111 0" 

"000001111111 0" 
"000011111111 0" 



9 BITS "0001111111110" 



10 Brrs "0011111111110" 



"0111111111110" 
"1111111111110" 



4 



7 
8 

9 

10 

11 

12 



WO 99/16224 



PCT/SE98/01598 



7/7 



Q Bit Allocation ^ 100 



Calculate for each subchannel a series 
of threshold values X correlated to a preset 
number of bits and store in modulator. 



T 



Fig. 10 



•102 



Obtain noise deviation estimate 0 and correct. 



•104 



Set number of active subchannels N S c to all subchannels. 



-106 



Compare 0|< to thresholds established for each subchannel and 
allocate corresponding number of bits \ for exceeded threshold 
to subchannel. For any subchannel where \ - 0, N sc = N sc -1. 




• 108 



Channel is insufficient for the 
targeted BER/Bitrate. 
System decides on step-down 
in BER or Bitrate demand based 
on service requirements. 



bits in 
than or equal 



greater 



number i 
excess of B tar g, 

he number of active- 
subchannels N sc ? 



YES 



116 



Subtract one bit form the 
bitallocation on all active 
subchannels. Any subchannel 
reduced to zero is inactive 
i.e. N sc = N sc -1 . 



Calculate distance d^ for each subchannel k. 



i 



• 118 



In order of decreasing distances, subtract one bit 
from each active subchannel until the sum of 
allocated bits is B target 



-120 




INTERNATIONAL SEARCH REPORT 



A. CLASSIFICATION OF SUBJECT MATTER 

IPC 6 H04L27/26 



Inter >nal Application No 

PCT/SE 98/01598 



According to international Patent Classification (IPC) or to both national classification and IPC 



B. FIELDS SEARCHED 



Minimum documentation searched (classification system followed by classification symbols) 

IPC 6 H04L 



Documentation searched other than minimum documentation to the extent that such documents are included in the fields searched 



Electronic data base consumed during the international search (name of data base and, where practical, search terms used) 



C. DOCUMENTS CONSIDERED TO BE RELEVANT 



Category 3 



Citation of document, with indication, where appropriate, of the relevant passages * 



Relevant to claim No, 



EP 0 753 947 A (ALCATEL BELL) 
15 January 1997 

see column 1, line 53 - column 2, line 4 
see column 2, line 54 - line 59 
see column 5, line 25 - line 56 

EP 0 753 948 A (ALCATEL BELL) 
15 January 1997 
see claim 1 

US 5 479 447 A (CHOW, CIOFFI) 
26 December 1995 

see column 3, line 57 - column 4, line 20 
see column 5, line 17 - line 23 
see column 7, line 11 - line 20 

-I— 



1-33 



1-33 



1-33 



m 



Further documents are listed In the continuation of box C. 



0 



Patent family members are listed In annex. 



" Special categories of cited documents : 

-A" document defining the general state of the art which Is not 

considered to be of particular relevance 
"E" earlier document but published on or after the International 

filing date 

V document which may throw doubts on priority claim(s)or 
which is cited to establish the publication date of another 
citation or other special reason (as specified) 

-0" document referring to an oral disclosure, use, exhibition or 
othermeans 

■p" document published prior to the international filing date but 
later than the priority date claimed 



T" later document published after the international tiling date 
or priority date and not In conflict with the application but 
cited to understand the principle or theory underlying the 
invention 

"X" document of particular relevance; the claimed invention 
cannot be considered novel or cannot be considered to 
involve an inventive step when the document is taken alone 

"Y" document of particular relevance; the claimed Invention 

cannot be considered to involve an inventive step when the 
document Is combined with one or more other such docu- 
ments, such combination being obvious to a person skilled 
In the art. 

document member of the aame patent family 



Date of the actual completion of the International search 



26 January 1999 



Date of mailing of the international search report 



02/02/1999 



Name and mailing address of the ISA 

European Patent Office, P.B. 5818 Patentlaan 2 
NL -2280 HVRiJswijk 
Tel (+31-70) 340-2040, Tx. 31 651 epo nl, 
Fax: (+31-70) 340-3016 



Authorized officer 



ScHven, P 



Form PCTYISA/210 (second sheet) (July 1 992) 



page 1 of 2 



INTERNATIONAL SEARCH REPORT 



Intel anal Application No 

PCT/SE 98/01598 



C.(Continuatlon) DOCUMENTS CONSIDERED TO BE RELEVANT 



Category a 



Citation of document, with indication, where appropriate, of the relevant passages 



Relevant to claim No. 



GB 2 300 546 A (SAMSUNG) 6 November 1996 
see page 5, line 21 - page 6, line 8 
see page 8, line 23 - line 26 



1-33 



Forni PCT/IS/V210 (continuation of second sheet) (July 1992) 



page 2 of 2 



INTERNATIONAL SEARCH REPORT 

.rformatlon on patent family members 



Inte onai Application No 

PCT/SE 98/01598 



Patent document 
cited in search report 



Publication 
date 



Patent family 
member(s) 



Publication 
date 



EP 075394/ 


A 

A 




.1 0G7 


All 

All 


DootbyD A 


Ol A1 1 Aft"? 

Z3-Q1-1997 










CA 


2180973 A 


12-01-1997 










US 


5812599 A 


22-09-1998 


EP 0753948 


A 


15-01- 


•1997 


AU 


5832896 A ; 


23-01-1997 










CA 


2180972 A 


12-01-1997 










US 


5790550 A 


04-08-1998 


US 5479447 


A; 


26-12- 


•1995 


NONE 






GB 2300546 


A 


06-11- 


•1996 


JP 


8307385 A 


22-11-1996 










US 


5832387 A 


03-11-1998 



Forni PCT/lSA/2ip (patent tamly annex) (July 1992) 



